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We use the recently proposed Kaluza-Klein (KK) reduction over the time dimension, within an 
effective field theory (EFT) approach, to calculate the next to leading order (NLO) gravitational 
spinl-spin2 interaction between two spinning compact objects. It is shown here that to NLO in 
the spinl-spin2 interaction, the reduced KK action within the stationary approximation is sufficient 
to describe the gravitational interaction, and that it simplifies calculation substantially. We also 
find here that the gravito-magnetic vector field defined within the KK decomposition of the metric 
mostly dominates the mediation of the interaction. Our results coincide with those calculated in 
the ADM Hamiltonian formalism, and we provide another explanation for the discrepancy with the 
result previously derived within the EFT approach, thus demonstrating clearly the equivalence of 
the ADM Hamiltonian formalism and the EFT action approach. 
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I. INTRODUCTION 



A novel effective field theory (EFT) approach treating the post-Newtonian (PN) formalism of General Relativity 
(GR) was introduced recently by Goldberger and Rothstein [Ij. It is very advantageous in applying the efficient 
standard tools of Quantum Field Theory to GR, notably handling the regularization required for higher order correc- 
tions in the PN approximation with the standard renormalization scheme. Moreover, an EFT approach enables us 
to handle a variety of physical situations, which are characterized by having multiple typical length scales. Initially, 
it was implemented on the evolution of an inspiralling binary to yield predictions of gravitational radiation [ll, . 
Later, it was used to obtain higher order thermodynamic results for higher dimensional Kaluza-Klein (KK) black 
holes fsj. Recently, it was further improved, and used to obtain first thermodynamic properties of higher dimensional 
KK rotating black holes (4)] . 

In this paper, we make a first application in the PN approximation of the KK reduction over the time dimension, 
proposed recently in ^ within the EFT approach, to obtain the next to leading order (NLO) spinl-spin2 interaction 
of a binary of spinning compact objects. The KK reduced action in the stationary approximation is shown here to 
be sufficient to describe the gravitational interaction to NLO in the spinl-spin2 sector. It is also shown to greatly 
simplify calculation, as well as provide further physical insight on the mediation of the interaction with spin. Our 
results coincide with those calculated in the ADM Hamiltonian formalism Q, demonstrating clearly its equivalence 
with the EFT action approach. Here, we provide another explanation for the discrepancy with the result previously 
derived within the EFT approach in Q, which with the addition of Q, was shown in the later to be canonically 
related to the ADM result. 

After the completion of this work, Q appeared which partially overlaps it. 



II. SPINNING OBJECTS IN AN EFT APPROACH WITH KALUZA-KLEIN REDUCTION 



To the order that we are calculating here, the action describing the dynamics of two spinning objects is given by 
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where Ta and Aq are the proper time and worldline parameters of the a'th particle's worldline jioj . The first term is 
the relativistic point particles (PP) interacting with gravity, where the PP spin coupling appears as the second term 
here, fi^i/ is the particle generalized angular velocity, expanded around flat spacetime in terms of tetrads and the 
connection of the metric, with the spin variable S^'^ conjugate to it. The last part of the action is the gravitational 
field interaction, including the usual Einstein-Hilbert (EH) action together with a gauge fixing term. Other higher 
order terms in the action do not contribute to the calculation of the NLO spinl-spin2 interaction as explained in 
The conventions used here are c = 1, a (-I-, — , — , — ) signature, and R^^^^p ~ ^^^^p ~ ^P^ua + ' ' '■ 

In the expansion of the metric around flat spacetime the metric is decomposed into potential and radiation modes, 
as explained in Both modes have the same typical time variation scale set by the binary orbital frequency v/r, 
but there is a different typical length scale for each, r and r/v respectively, where r and v are the typical orbital 
parameters of the binary, and we are working in the PN limit v ^ I. Thus, as far as the potential gravitons are 
concerned, the radiation component is just a slowly varying background held of soft momenta gravitons. Moreover, 
considering these typical scales, we see that the potential gravitons are off shell, with their frequency being much 
smaller than their momentum, and thus can be approximated as stationary to leading order. 

Based on these observations a KK reduction over the time dimension of the potential field was suggested and 
used to simplify the gravitational action, and consequently the EFT calculation 0,0] ■ At first stage, the metric is 
parametrized according to the Kaluza-Klein ansatz 
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This just defines a set of new fields (</>, ^^,7^), i,j — 1,2,3, the scalar field <j), which corresponds to the Newtonian 
potential, the gravito-magnetic vector field Ai, and the 3 dimensional symmetric 2-tensor field 7^^, the non-relativistic 
gravitational (NRG) fields as discussed in Q. Next, we suppress the time dependence of the fields to obtain the KK 
reduced action for the gravitational field, given by 



Skk — 



(3) 



where -F,-, 



d^A, 



djA, 



After obtaining the KK reduced action the gauge fixing term is to be chosen to set 



completely the gravitational action. 
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Thus, it remains to define the different NRG graviton propagators, which are used in our calculations. The 
propagators can be read off from the terms that are quadratic in the fields in the gravitational action, depending on 
the choice of the gauge fixing term of course. We use here the natural gauge adequate for the NRG fields, namely the 
Lorentz gauge for the vector field Ai, and the linear harmonic gauge for the 2-tensor field 7^^- in 3D. Thus, the gauge 
fixing term reads 



(4) 



Quadratic terms with time derivatives, such as that we shall see below in Eq. ([M)) . are suppressed as subleading 
corrections in powers of v'^ as explained above: the potential graviton propagators are instantaneous within the 
leading stationary approximation, representing off shell gravitons. The NRG scalar, vector, and tensor graviton 
propagators in the harmonic gauge are thus given by 

(7.,(xi)7fe/(x2)> = /iC^T^e^"^-^"^-"^^ P^rM , 

where Pij-ki = \ {Sik^ji + SuSjk — 25ijSki). Note that the propagator of the gravito-magnetic vector field, which is the 
most common in this interaction as we shall see, is diagonal in the harmonic gauge - an advantageous feature of the 
NRG parametrization. 

Before we specify the calculation, we need to consider the gauge for the redundant spin degrees of freedom. We 
may choose the Newton- Wigner (NW) coordinates, which yield the spin supplementary condition (SSC) given by [12] 

mS'O'^ = S^"'p^. (6) 

We may also choose a conjugate gauge for the body fixed tetrads given by 

so that in the rest frame we reduce to pure rotations as physically required [l^]. Using these constraints we get 

5-^0 ^ ^g^J^J ^ ^S'^Aj + O(w^), (8) 

where the subleading temporal spin entries S'**' in the local Lorentz frame depend both on the coordinate velocity 
and on the metric field. Note the 0{v^) subleading term appearing here, which is absent in flat spacetime, where 
the local Lorentz and coordinate velocities coincide, and only terms explicitly odd in the power of the velocity exist. 
Note also that it is the gravito-magnetic vector field, which is involved in this leading curved spacetime contribution 
to the temporal spin entries. 

We may also choose the covariant SSC, given by S^'^Pi, — [T^]- However, to the order considered here this is 
equivalent to the condition given by S^'^Ui, = 0, i.e. in the covariant SSC, we have for the spin temporal components 

S^a ^ g^J^J ^ ^S'JAj +■■■, (9) 

which again implies a suppression of the temporal spin components 5*" by one order of v with respect to the spatial 
components S^K However, here there is no infinite expansion explicit in v, yielding contributions to all orders, as 
opposed to Eq. dS]). These constraints, e.g. those expressed in Eqs. ((SUH]), eliminate the redundant unphysical degrees 
of freedom, so that the spin may be represented by a 3-vector defined as S"^^ = e^^^S^ . 

We go on to expand the spin part of the action in terms of the NRG fields to obtain the following leading order 
(LO) spin graviton coupling 

L2 = \S''F,„ (10) 



where the subscript is standard notation for the power of the orbital velocity of this term in the action. Thus, the LO 
spinl-spin2 potential follows from a one-graviton exchange of the gravito-magnetic vector graviton, with two of these 
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FIG. 1. Feynman diagram of the leading order spinl-spin2 interaction. The thick sohd hnes represent the time evolution 
of the point particles worldlines. The blobs represent spin insertions on the worldline. The dashed line represents the Ai 
gravito-magnetic vector field propagator. 



LO spin vertices, and will scale accordingly as (2PN). The corresponding Feynman diagram is depicted in Fig. 1, 
resulting in the weW knovifn LO spinl-spin2 potential given by 



= -^01 ■S2~3Si-nS2- n), 



(11) 



where r = jxi — X2I, and n = ^. 

For the next order we must also take into account the subleading temporal spin entries (i.e. 5'^}, where we are 
considering the SSC on the level of the action , which is accurate to the order we are calculating here . Considering 
vertices involving spin to order v^, we have 



(12) 



where the leading order term of the temporal spin entry 5°' is considered in the second term here. Calculating to 
order (3PN), we should include diagrams with two insertions of as depicted in Figs. 2(al), 2(a2). The values of 
these diagrams are given by 



Fig. 2(al) 
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Fig. 2(a2) 



G r 



Si ■ S2 vi ■ V2 — Si ■ V2 S2 ■ vi — 3Si X vi ■ n S2 ^ V2 ■ n 

^'n'S°'n^ 

S^' ((5^2 X V2y - 3n' S2XV2-nj+ S°' ((5i x vi)' - 3n' Si x vi ■ n 
~2Si ■ S2 vi ■ V2 - Si ■ Hi S2 ■ V2 + 2Si ■ V2 82- vi+ iSi ■ n S2 ■ n vi ■ V2 
+6S1 X vi ■ n S2 X V2 ■ n — S^i x V2 ■ n S2 x vi ■ n 



(13) 



(14) 



Here, a multiplicative factor of J dt is suppressed in the values of diagrams, and is to be further omitted from all 
diagram values. 

For the NLO we should also include time derivatives, thereby departing from stationarity. At order we must 
include 



(15) 



Both terms contain temporal spin entries S^^ considered in their LO. Note the second term here, first containing time 
dependence of the potential field, that of the gravito-magnetic vector. Thus, we should also include diagrams with 
one insertion of and the LO spin vertex L2 as depicted in Fig. 2(b), which equals 



Fig. 2(6) = 



G 
G 



35?^ {n'S2 xvi-n- {S2 x n)Vi • nj + 3S^' (n'Si x V2 ■ n - {Si x n)V2 • nj 
dtS"i'iS2 X nf -dtSl\Si X n)' 



(16) 



However, note that due to the time derivative term this diagram may be evaluated in two ways. The one we took 
here is straightforward, but yields acceleration terms. Such terms may seem undesirable since they call for the use of 
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the equations of motion (EOM) at the level of the Lagrangian, a procedure which is known to be incorrect in many 
(in)famous examples. However, as originally noted by [14| . and later formally evolved and treated, e.g. in 15|, a 
substitution of low order EOM in higher order terms in the level of the Lagrangian is a correct procedure in GR, and 
is equivalent to performing a coordinate transformation. For example, if we consider the NW SSC, i.e. Eq. (jS]) at the 
level of the action, we evaluate Fig. 2(b) according to Eq. (IT6l) as 

G [ ^ - ■ - ^ -> 

-ZSi ■ S2 {{vi ■ nf + {v2 ■ nf) + 3(S'i • V2){S2 ■ n){v2 ■ n) + 3(5*2 • vi){Si ■ n){vi ■ n) 



Fig. 2(6) 



2r3 



+3(S'i X vi) ■ n{S2 X vi) ■ n + 3(51 x V2) ■ n{S2 x V2) ■ n 
G 



2r2 



G 
2^ 



{Si X ai) ■ {S2 X n) - (52 x ^2) • {Si x n) 
iSi ■ S2 {{vi ■ nf + {V2 ■ nf) + 3(5i • V2){S2 ■ n){v2 ■ n) + 3(^2 • vi){Si ■ n){vi ■ n) 

+3(S'i X vi) ■ n{S2 X vi) ■ n + 3(5i x V2) ■ n{S2 x V2) ■ n 
G2(toi + TO2) 



2r4 



Si ■ S2 — Si ■ n S2 ■ n 



Here, we eliminated the acceleration terms using the LO EOM, given by 

G{mi + 7712) ^ 



a = ai — a2 = 



(17) 



(18) 



resulting in the nonlinear term in G appearing in the last line of Eq. (I17p . By this substitution we implicitly have 
coordinate transformations made (leading us to the result in ADM coordinates). 

Alternatively, the time derivative term from Eq. (jl5p may be evaluated by flipping the time derivative between the 
two PP worldlines, namely by using the identity 

dtidt2 dt,S{ti - t2)f{ti)g{t2) = - j dtidt2 dt25{ti - t2)f{ti)g{t2). (19) 

Such an evaluation of the diagram in Fig. 2(b) avoids the acceleration terms, and the use of EOM. Again, considering 
the NW SSC at the level of the action, it yields the following seemingly different value for the diagram 



Fig. 2{bU. = ^ 



25i ■ S2 vi ■ V2 ~ 2Si ■ V2 S2 ■ vi + Si ■ vi S2 ■ vi + Si ■ V2 S2 ■ V2 



(al) 



(a2) 



v'(\ Qv' v'Ch -(g)- -Ov 



(b) 



(c) 



FIG. 2. Next to leading order spinl-spin2 interaction Feynman diagrams of one-graviton exchange. The solid line represents the 
scalar field propagator. The double line represents the 2-tensor field propagator. The cross vertex corresponds to an insertion 
of the graviton kinetic term. Diagram (b) should be included together with its mirror image. 
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-5^1 • 5^2 {{vif + {V2)^) - 6(5i • S2){vi ■ n){v2 ■ n) 
-3{Si ■ V2){S2 ■ n){vi ■ n) + 3(52 ■ vi){Si ■ n){v2 ■ n) 
|-3(S'i X vi) ■ n{S2 X vi) ■ n + 3(5*1 x W2) • n{S2 x V2) ■ n 



(20) 



This was the route taken in ?! , thus requiring an explicit canonical transformation to reproduce the equivalent result 
in the ADM coordinates after an addition of a missing piece in the potential as shown later in [s^ . The time derivative 
term of Eq. p5|) also yields spin precession terms, thus similarly requiring the use of the LO spin EOM, which scales 



as 



m 



dS_ 
'dt 



Gm 



Sv. 



(21) 



This yields a 4PN contribution, and can therefore be taken here as zero. 

There is a further contribution from the one-graviton exchange sector, which arises from the leading correction to 
the potential graviton propagators, i.e. the graviton kinetic terms in the gravitational action. It is a first departure 
from the stationarity approximation in the gravitational action. Extracting from the EH action the quadratic time 
dependent terms, which are relevant to the order we are calculating here, we have 



39' ^ L_ 



dtcPx 



1 



1 



(22) 



Thus, in order to reduce the number of contributing diagrams as much as possible, we refine the gauge fixing term of 
Eq. Q to eliminate undesired quadratic vertices. Hence, we reset it to be 



Sgf — 



1 



327rG 



so that to the order considered here, we are left with the following time dependent quadratic term 

1 



Seh + Sgf ^ — 
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d x{doAi 



(23) 



(24) 



Hence, the leading deviation from stationarity in the gravitational self interaction in the spinl-spin2 sector involves 
again the gravito-magnetic vector field alone. However, note that since the gravito-magnetic vector has no kinetic 
term in the original EH action, this contribution arises as a pure gauge fixing correction, i.e. it does not involve 
corrections to the KK reduced action. This last contribution is depicted in the diagram shown in Fig. 2(c), which 
equals 



P^9. 2(c) = ^ 



-5*1 ■ S2 Vi ■ V2 + Si ■ Vi S2 ■ V2 + Si ■ V2 S2 ■ Vi 

+3{Si ■ S2)ivi ■ n){v2 ■ n) - 3(51 • n){S2 ■ n){vi ■ V2) 
-3(5i • vi){S2 ■ n){v2 ■ n) - 3(^2 • V2){Si ■ n){vi ■ n) 
-3(5i • V2){S2 ■ n){vi ■ n) - 3(^2 • vi){Si ■ n){v2 ■ n) 
+ 15(5i • n){S2 ■ n){vi ■ n){v2 ■ n) . 



(25) 



Here, one also makes use of the LO spin EOM to eliminate precession terms. 

At order 3PN nonlinear contributions in G must also be included. First, we have the two-graviton exchange 
(including only connected diagrams after stripping off worldlines). The relevant terms in the Lagrangian, which are 
quadratic in the metric, are given by 



T3 

^4 



S^'F.,, 



\s''A,d,^ 



(26) 



However, it is important to note that this sector contains a further contribution from the subleading 0{v'^) term 
of the 5"* spin entries, e.g. in Eqs. ([8l[9|), which incorporates the effect of curved spacetime. It arises as the NLO 
substitution in the 5°' term in Eq. p^ . For example, if we consider the NW SSC, it equals 



T3 

^4:(SSC) 



1 



(27) 
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(b) 



FIG. 3. Feynman diagrams of nonlinear spinl-spin2 interaction in next to leading order. The black blobs represent mass 
vertices on the worldline. These diagrams should be included together with their mirror images. 



This contribution was considered only in [sj, and then regarded as an additional spin-orbit effect, but here it is 
incorporated into the Feynman diagrams as being a proper spinl-spin2 contribution. Hereby, we consistently consider 
the SSC at the level of the action, a procedure accurate to the order calculated here as we noted already, while not 
being generally valid. The terms in Eqs. (|26[ \Tll scale as w*, so that diagrams containing them must also include a 
LO spin insertion, and a LO mass insertion from the non-spinning part of the PP action given by 



Lq = —m(j). 

This diagram is shown in Fig. 3(a), and considering Eq. (|26p . it equals 

G^(toi + TO2) 



Fig. 3(a) 



55*1 • 5*2 — 13Si ■ n S2 ■ ft 



In addition, the contribution from e.g. Eq. (I?7l) . arising from the NW SSC, yields 

G^(mi + 7712) 



Fig. i{a)ssc 



2r4 



Si ■ S2 — Si ■ n S2 ■ n 



(28) 



(29) 



(30) 



Finally, we must include the contribution from those diagrams including the three-graviton vertex, which scales as 
To the order considered here diagrams including this vertex should therefore include two LO spin insertions 
and one LO mass insertion, with the 3-graviton vertex given by the following cubic term from the KK reduced action 

1 



Skk ^ 



SttG 



dtd-^x [(j>^^AJ {d^Aj - djA,)] . 



(31) 



This vertex is easily read from the KK reduced action, already making computation much faster compared to the 
tedious extraction of the cubic part of the EH action, which contains a very large number of terms, each with a 
complicated tensor index structure as in The corresponding diagram is depicted in Fig. 3(b). This diagram 
contains a loop integral, and requires regularization. It is handled with dimensional regularization by the usual 
techniques, see e.g. [16j. The diagram is evaluated as 



Fig. 3(6) 



2G2(toi +7712) 



Si ■ S2 ~ 25*1 ■ fi S2 ■ n 



(32) 



The long calculation required in Q, calling for the aid of a symbolic manipulation software, becomes extremely feasible 
by hand as noted already in Q. Both diagrams of the nonlinear interaction are included with their mirror images. 



III. NEXT TO LEADING ORDER SPIN1-SPIN2 POTENTIAL 



Summing up all of the contributions from the diagrams, we obtain the NLO spinl-spin2 Lagrangian. Here, we give 
the result without substituting in the SSC, leaving the 5'°' entries as independent degrees of freedom, and without 
taking into account the SSC field contribution e.g. in Eq. (I?7l) . Next, we will consider both the implementation 
of NW SSC at the action level, and the covariant SSC at the Hamiltonian level. The results are to be compared 
with Eq. (2.11) of [6J, and we will show that both ways lead to an agreement with the result obtained in the ADM 
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Hamiltonian formalism. Our Lagrangian then reads 
G 



-NLO 
'S1S2 



2r3 



5*1 • 5*2 (3wi • V2 ~ 3vi ■ nv2 ■ n) + Si ■ V1S2 ■ ^2 - 3Si ■ V2S2 ■ vi 
—3Si ■ nS2 ■ n {vi ■ V2 + 5vi ■ nv2 ■ n) + iSi ■ V1S2 ■ nv2 ■ n + 3S'2 • V2S1 ■ nvi ■ ft 
+3S'i • V2S2 ■ nvi ■ n + 'SS2 ■ viSi ■ nv2 ■ n — 6(S'i x vi) ■ n{S2 x V2) ■ n + 6{Si x V2) ■ n{S2 x vi) ■ ft 
G^(toi + 7712) 



-3 



Si ■ S2 ~ 35*1 • nS2 ■ n 



G r 



G 



S^^'Sl' - iSl'n'Slhi^ + Si' {{S2 X V2Y - 3n'S2 x V2 ■ nj + S^' (^{Si x vi)' - Sn'Si x vi ■ ftj 
3Sl' (n'S2 xvi-n~ {S2 x n) Vi • nj + 3S^' (n'Si x V2 ■ n - {Si x n)V2 • nj 

dtSl\S2 X ny - dtSfiSi X ny] , 



(33) 



where the SSC dependent part appears in the last three lines of the above expression. 



A. Implementation of NW SSC at the Action level 

If we consider the NW SSC at the level of the action, i.e. we substitute in the 5*" entries from Eq. ([5]), and take 
into account the contribution from Eqs. ((27l [30|) . we obtain 



j^NLO 

S1S2 



G 
'2^ 



Si ■ S2 (^^^1 • ^2 - 3^1 • nv2 -n + S {{vi ■ nf + {v2 ■ nf)^ + 5^1 • V1S2 ■ V2 - • V2S2 ■ vi 

— 35*1 • nS2 ■ n (vi ■ V2 + 5vi ■ nv2 ■ n) + 35i • V1S2 ■ nv2 ■ n + 3S2 ■ V2S1 ■ nvi ■ n 

+iSi ■ V2S2 ■ nvi ■ n + 3S2 ■ viSi ■ nv2 ■ n — 3Si ■ V2S2 ■ nv2 ■ n — 3S'2 • viSi ■ nvi ■ n 

3 ^ 

— 3(S'i X vi) ■ n{S2 X vi) ■ n — 3(S'i x V2) ■ n{S2 x V2) ■ n + 2 ("^i ^ '^i) ' ""("^z x V2) ■ n 



+&{Si X V2) ■ n{S2 X vi) -n 
G^(mi + m2) 



G 
2r2 



{Si X di) ■ {S2 X n) - {S2 X 0,2) ■ {Si X n) 



2r3' 



2r4 

G'^{mi + TO2) 



5Si ■ S2 ~ 175i • nS2 ■ n 



Si ■ S2 ~ 3Si ■ nS2 ■ n 



(34) 



where in the last equality we have used the LO EOM for the accelerations. 

We should make a Legendre transformation with respect to the coordinate velocities in order to obtain the NLO 
spinl-spin2 Hamiltonian. The relevant part of the Lagrangian reads 



LO 



~ yso 

where is the Newtonian Lagrangian given by 



LO 



S1S2 



NLO 
'S1S2 ' 



2 

1 ^2 , Gmim2 
Ln = x2^maw, + , 



(35) 



(36) 



VgQ is the LO spin-orbit potential, given by (e.g. from 6]) 



yLO 

^so 



Gm2 



Si-{vixn~ 2v2 xn) + Sfn' +[1^2] 



Gm2 



Si 



vi X n - 2v2 X n] + [I ^ 2], (37) 



where in the second equality the NW SSC from Eq. ([8]) is implemented. Vg^^ is the LO spinl-spin2 coupling given 
here in Eq. (|lip . and L^^s^ is the NLO spinl-spin2 coupling that we calculate here. Note that the LO spin-orbit 
potential already requires the use of SSC, so we need to pay special attention to it. The canonical momenta are given 
by 



Pa 



dL 

dVa 



(38) 
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so that e.g. 



G /3 



pi = mivi - ^ 



Consequently, we have for the velocity, e.g. 



. Pi ^ G 

vi = \ 5 

mi mir^ 



-TO2(S'i X r?) + 2mi(S'2 x n) + 



-TO2(S'i X n) + 2toi(S'2 X n) 



5^ 



ml 



(39) 



(40) 



The Hamiltonian is obtained by substituting in the velocities given in terms of canonical momenta in the following 
expression 



H = ^ Va-Pa- L. 



(41) 



Note that the 



terms get canceled, and do not contribute here. 



The Hamiltonian then takes the following form 

H = Ho + Hso + Hss, 



(42) 



where Hq is a part of the orbital Hamiltonian, Hso is a part of the spin-orbit (linear in spins) Hamiltonian, and Hss 
is the spin-spin (quadratic in spins) sector of the Hamiltonian up to NLO. Thus, the NLO spinl-spin2 Hamiltonian 
equals 



H 



NLO _ tNLO , 3G'2(mi+m2) 



S1S2 



Si ■ S2 ~ Si ■ nS2 ■ n 



(43) 



and finally we obtain 



H 



NLO 
S1S2 



G 
2^ 



— ^ — ( \Si ■ S2P1 ■ P2 - 35^1 • S'2Pi ■ np2 ■ n + Si ■ P1S2 ■ P2 - 7:^1 ■ P2S2 ■ Pi 
TO1TO2 \2 2 

—35*1 • nS2 ■ npi ■ P2 — 155i • nS2 ■ npi ■ np2 ■ ft + 35i • P1S2 ■ np2 ■ n + 3S2 ■ P2S1 ■ npi ■ ft 

3^ ^ ^ ^ 

+35*1 •_P2'5'2 • fipi ■ ft + 35*2 • piSi ■ ftp2 ■ fi+ - {Si X pi) • fi{S2 x P2) ■ ft + 6{Si x P2) ■ ft{S2 x pi) ■ fi 



3 
m 



2 [Si ■ S2{pi ■ flY - S2 ■ PiSi 
1 ^ 



ftpi ■ ft — [Si X pi) • fi{S2 X Pi) - ft 



H — o (Si ■ S2{p2 ■ fif - Si ■ P2S2 ■ ftp2 ■ fi - 01 X P2) ■ fl{S2 X P2) ■ ft 
mi V 



G2(mi + m2) 



6Si ■ S2 — 125i • fiS2 ■ fi 



(44) 



This is just equal to the result given in Eq. (2.11) of [6*], obtained by the ADM Hamiltonian formalism. 



B. Implementation of the Covariant SSC at the Hamiltonian level 

Let us consider our NLO spinl-spin2 interaction Lagrangian given in Eq. p3p with the implementation of the 
covariant SSC only at the Hamiltonian level, i.e. considering the 5*° as independent degrees of freedom until the 
obtainment of a Hamiltonian, and only then implementing the covariant SSC. 

Again, we need to Legendre transform Eq. ([55)1 . but now instead of Eqs. ([5^ HO)) , we will have from Eq. ([57]) e.g. 

pi = mivi - {m2{Si xfi) + 2mi{S2 x n)j + , (45) 

vi^ — + ^ (m2(5i X n) + 2mi(52 x n)) - (46) 

mi mir"^ V / ml 
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Thus, Eq. p3l) is Legendre transformed into the fohowing NLO spinl-spin2 Hamihonian 



rrNLO_ tNLO , 2G'^(toi +TO2) 



Si ■ S2 ~ Si ■ nS2 ■ n 



G 



— [rn^Sfv} + 2)] 



G 
2^ 



1 



Si ■ S2 (3pi • p2 - 3pi ■ np2 ■ n) + Si ■ P1S2 ■ P2 - 3Si ■ P2S2 ■ Pi 



_miTO2 

-35^1 • nS2 ■ n{pi ■p2 + 5pi ■ np2 ■ n) + 3Si ■ P1S2 ■ np2 ■ n + 3^2 • P2S1 ■ npi ■ n 
+3Si ■ P2S2 ■ npi ■ n + 3S2 ■ PiSi ■ np2 ■ n - 6(S'i x pi) ■ n{S2 x P2) ■ n + 6{Si x P2) ■ n{S2 x pi) ■ n 



G 



SfSl' - iSfTlSl^n^ + — ( (S-a X ^2)' - 'in'-S2 x p2 ■ n] ^ [{Si x pif - 3n*5'i xpi-n] 
m2 \ / mi V / 



Oi 



+ 

+ 



yi^S2 X pi ■ ft — {S2 X HYpi ■ nj H yn^Si x p2 ■ n ~ (Si x r?)*p2 • n 



3Sl 
mi 



55*1 • 5*2 — llSi ■ nS2 ■ n 



G 



G 



dtS"i'iS2 X nf - dtS^'iSi X rty + — [m2Sf'n' + (1 o 2)] 



(47) 



Note the last term here originating from the SSC dependent part of the LO spin-orbit interaction in Eq. ([37)) . We 
need to keep track of this term, as it contributes to the NLO spinl-spin2 interaction too. Next, we go on to insert the 
covariant SSC to our Hamiltonian, using Eqs. ([9]), (|46p . and the metric for spinning binary black holes in harmonic 
coordinates, e.g. as in 17[. It is well known that the variables then obtained will be non-canonical. However, we shall 
proceed to make non-canonical transformations that relate them to the canonical variables . 
Hence, upon the insertion of the covariant SSC our non-canonical Hamiltonian reads 



H 



NLO_ 
S1S2 - 



G 
2r3 



1 



miTO2 



-Si ■ S2 (3pi • p2 + 3pi ■ np2 ■ n) + Si ■ P1S2 ■ P2 + 35i • p2<5'2 • Pi 



-35*1 • nS2 ■ ri{pi ■p2 + 5pi ■ np2 ■ n) 

-1-35*1 ■ Pi<S'2 • np2 ■ fi + 3S2 ■ P2S1 ■ fipi ■ n + 3Si ■ P2S2 ■ npi ■ n + 35*2 • piSi ■ np2 ■ n 
-|-12(5i X pi) ■ n{S2 X P2) ■n + 6{Si x P2) ■ n{S2 x pi) ■ nj 
6 



H 2 (^i ■ '^2(pi ■ nf - 52 • _pi5i ■ npi-n - {Si x pi) ■ n{S2 x pi) ■ ft 

+ ~~2 (-^i ■ ^2{P2 ■ nf - Si ■ P2S2 ■np2-n- {Si x pa) • n{S2 x P2) ■ 

77T2 ^ ' 

G'^{mi + 1112) 



7Si ■ S2 — 135i • n52 • n 



(48) 



We note that had we implemented the covariant SSC at the level of the action, we would have arrived at the same 
Hamiltonian appearing here. Now, we go on to apply the following non-canonical transformations of variables 



5i 5i 

r 1 — !• fi + 



Pi ^ Pi ■ Si ^ _ ^ Pi X {Si X pi) 



2m{ 2mi 

Si X pi G 
2ml 



2mf 



mir^ 



S2 X n] X Si 



(49) 



(50) 



These transformations are just the generalization of the mapping between the covariant and NW SSC variables in flat 
spacetime [l^l- The transformation of the spin variable in (H^ is just similar to that of flat spacetime, whereas for 
the center of mass coordinate in (|50p . a higher order PN transformation, generalized for curved spacetime, is required 
[isj . Finally, after the application of these transformations, we arrive at the canonical Hamiltonian of 6] obtained by 
the ADM Hamiltonian formalism, appearing here in Eq. (HH). 



C. Summary of Results 

Our results coincide with those of @ obtained by the ADM Hamiltonian formalism, and are also equivalent to those 
previously calculated in the EFT approach in , after taking into account the addition appearing in [8| . In the later 
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Q, the equivalence to the result in [sj was shown using a canonical transformation relating the results. However, the 
origin of the discrepancy between the results was not completely clear before. First, there is the addition appearing 
in Our results here in Sec. IIII Al confirm that to NLO in the spinl-spin2 interaction applying consistently the SSC 
at the level of the action in the EFT approach is a valid procedure. Here, the additional term from Eq. ([5(7|) . regarded 
in ^] as a spin-orbit effect, is incorporated in the Feynman diagrams, and is thus included as a proper spinl-spin2 
Feynman diagram - that of Fig. 3(a). By that we consistently applied the SSC at the level of the action. While in 
general applying the SSC at the action level is not a valid procedure, it is correct to NLO in the spinl-spin2 interaction 
due to power counting considerations. 

Second, we noted here in the evaluation of the diagram in Fig. 2(b), where the time derivative term in the spin 
graviton coupling of Eq. (IT5t is present, that there are two ways to evaluate the diagram. The one taken here in 
Eqs. (IT6l [TTt is straightforward, but yields acceleration terms in the Lagrangian. Though one may be reluctant to 
have such terms, calling for the use of EOM in the Lagrangian, which is known to be incorrect in many examples, 
the substitution of low order EOM into higher order terms in the level of the Lagrangian is known to be a correct 
procedure in GR, and it is equivalent to a coordinate transformation. Indeed, we eliminated the acceleration terms 
by substitution of the LO EOM, and obtained the nonlinear contribution in G in Eq. (|17p . which appeared to be out 
of place in the one-graviton exchange sector. Hereby, we implicitly had a coordinate transformation made. However, 
as we noted, there is a second possible way to evaluate this diagram, which avoids acceleration terms, by flipping 
the time derivative involved between the two particles worldlines. This was the route taken in [7], thus requiring the 
canonical transformation in [Sj], in order to get the result of @ in the ADM coordinates. 

To conclude, our results demonstrate exactly that the position, velocity, and spin variables used in the EFT 
approach relate to the canonical ones via standard coordinate transformations and variable redefinitions, as opposed 
to the reservations expressed in |6i]. 

IV. CONCLUSIONS 

In this paper we made a first application in the PN approximation of the KK reduction over the time dimension, 
proposed in Q within the EFT approach, to calculate the next to leading order gravitational spinl-spin2 coupling 
of a binary of spinning compact objects. The KK decomposition applied to the metric to yield the NRG fields, 
and the KK reduction applied to the EH action, are shown here to improve both the calculation and the physical 
interpretation of the interaction substantially. On the more technical level, all vertices simplify, most notably the 
previously voluminous and complicated 3-graviton vertex is exchanged with the simple (jjF"^ cubic vertex of the KK 
reduced action. Moreover, the propagator of the gravito-magnetic vector, which is the most common in this spin 
interaction, is diagonal, thus facilitating most calculations, and the handling of polarization. On the physical level, 
the reduced KK action is shown here to be sufficient for the description of the gravitational self interaction within 
the stationary approximation to NLO in the spinl-spin2 interaction, with no need for corrections to it. We also find 
here that similarly to the LO spinl-spin2 interaction mediated exclusively by the gravito-magnetic vector graviton, 
the NLO spinl-spin2 interaction is also mostly mediated by the vector graviton, with the two first deviations from 
stationarity, and the leading metric correction term in the spin temporal entries associated with this vector graviton. 
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